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Abstract 

A model of a strongly correlated electron liquid based on the fermion condensation (FC) 
is extended to high-temperature superconductors. Within our model, the appearance of FC 
presents a boundary separating the region of a strongly interacting electron liquid from the 
region of a strongly correlated electron liquid. We study the superconductivity of a strongly 
correlated liquid and show that under certain conditions, the superconductivity vanishes at 
temperatures T > T c ~ T no d e , with the superconducting gap being smoothly transformed 
into a pseudogap. As the result, the pseudogap occupies only a part of the Fermi surface. The 
gapped area shrinks with increasing the temperature and vanishes at T = T*. The single- 
particle excitation width is also studied. The quasiparticle dispersion in systems with FC 
can be represented by two straight lines characterized by the respective effective masses Mt c 
and M£, and intersecting near the binding energy that is of the order of the superconducting 
gap. It is argued that this strong change of the quasiparticle dispersion at the binding can be 
enhanced in underdoped samples because of strengthening the FC influence. The FC phase 
transition in the presence of the superconductivity is examined, and it is shown that this 
phase transition can be considered as kinetic energy driven. 

PACS numbers: 71.27.+a, 74.20.Fg, 74.25.Jb 

1 Introduction 

Unusual properties of the normal state of high-temperature superconductors have been attracting 
attention for a long time. In describing these properties, which are well beyond the standard Fermi 
liquid theory, the notion of a strongly correlated liquid has emerged (see e.g. 0, |[). Later on, 
angle-resolved photoemission studies revealed unusual properties observed in underdoped samples, 
with the leading edge gap discovered up to the temperature T* > T c . This behavior is interpreted 
as coming from the pseudogap formation; it was observed in a number of underdoped compounds 
such as YBa2Cus06+x, Bi 2 Sr 2 CaCu208+5, etc. As T increases above T*, the pseudogap closes, 
leading to a large Fermi surface and an extremely flat dispersion in electronic spectra, which is 
called the extended Van Hove singularity [3 — 7] . Thereupon recent discovery of a new energy scale 
for quasiparticle dispersion in Eh^S^CaC^Os+i in the superconducting and in normal states f| [J 
have brought new insight to the physics of high temperature superconductors (HTS), imposing 
serious constraints upon the theory of HTS. The newly discovered additional energy scale manifests 
itself as a break in the quasiparticle dispersion near (50 — 70) meV, which results in a drastic change 
of the quasiparticle velocity [§, |], [K| . Such a behavior is definitely different from what one could 
expect in a normal Fermi liquid. 

A correlated liquid can be described in conventional terms, assuming that the correlated regime 
is related with the noninteracting Fermi gas by adiabatic continuity. This is done in the well-known 
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Landau theory of the normal Fermi liquid, but the question arising at this point is whether this 
is possible. Most likely, the answer is negative. To tackle the above-mentioned problems, we 
consider a model where a strongly correlated liquid is separated from the conventional Fermi 



liquid by a phase transition related to the onset of the FC |TT], |12|1 . The purpose of our paper is 
to show that without any adjustable parameters, a number of fundamental problems of strongly 
correlated systems are naturally explained within the model. The paper is organized as follows. In 
Sec. 2, we consider the general features of Fermi systems with the FC. In Sec. 3, we show that the 
pseudogap behavior can be understood within the standard BCS superconductivity mechanism 
provided the appearance of FC is taken into account. In Sec. 4, we analyze the condensation energy 
that is liberated when the system in question undergoes the superconducting phase transition 
superimposing on the FC phase transition. In Sec. 5, we describe the quasiparticle dispersion and 
lineshape. Finally, in Sec. 6, we summarize our main results. 



2 The main features of liquids with FC 



We first consider the key points of the FC theory. The FC is related to a new class of solutions of 
the Fermi liquid theory equation [|i~3|j 
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6n(p } T) 
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for the quasiparticle distribution function n(p, T) depending on the momentum p and the tem- 
perature T. Here F is the free energy, /i is the chemical potential, and s(p,T) = SE/Sn(p,T) 
is the quasiparticle energy, which is a functional of n(p, T) just like the energy E and the other 
thermodynamic functions. Equation (1) is usually represented as the Fermi-Dirac distribution 



n(p, T) — < 1 + exp 



( g (p,T)-/i) 
T 



(2) 



In a homogeneous matter and at T = 0, one obtains from Eq. (2) the standard solution rip(p,T = 
0) = 6(pf — p), with e(p ~ pf) — ^ = Pf{p — Pf)/MI, where pp is the Fermi momentum and M£ 
is the commonly used effective mass HIT 



1 1 de(p, T = 0) 
Ml p dp 



P=Pf- 



(3) 



It is assumed to be positive and finite at the Fermi momentum pp. This implies the T-dependent 
corrections to M£, the quasiparticle energy e(p), and the other quantities start with T 2 -terms. 

But this solution of Eq. (1) is not the only one possible. There exist "anomalous" solutions 
of Eq. (1) associated with the so-called fermion condensation [|ll], [14], Being continuous and 



satisfying the inequality < n(p) < 1 within some region in p, such a solution n(p) admits a finite 
limit for the logarithm in Eq. (1) as T — > 0, yielding 



e(p) 



5n{p) 
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Equation (4) is used in searching the minimum value of E as a functional of n(p) under the 
assumption that a strong rearrangement of the single-particle spectrum can occur. We see from 
Eq. (4) that the occupation numbers n(p) become variational parameters: the solution n(p) 
exists if the energy E is decreased by alteration of the occupation numbers. Thus, within the 
region pi < p < pf, the solution n{p) deviates from the Fermi step function n F {p) such that the 
energy e(p) stays constant, while n(p) coincides with n F (p) outside this region. As a result, the 
standard Kohn — Sham scheme for the single-particle equations is no longer valid beyond the FC 
phase transition point [[HJ. This behavior of systems with the FC is clearly different from what 
one expects from the well known local density calculations; therefore, these calculations are not 
applicable to systems with the FC. On the other hand, the quasiparticle formalism is applicable 
to this problem, because as we see in what follows, the damping of single-particle excitations is 
not large compared to their energy JTH] . It is also seen from Eq. (4) that a system with the FC 
has a well-defined Fermi surface. 

It follows from Eq. (1) that at low T, new solutions within the interval occupied by the fermion 
condensate have the spectrum e(p,T) that is linear in T [[15], |j~7|j , 

e (p, T) - /i(T) * iP ~l F)PF * T(l - 2n(p)) « T f . (5) 
Here Tf is the quasi-FC phase transition temperature above which FC effects become insignificant 
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?L„tizA„**<L (6) 
e F 2Me F Q F ' v ; 

where M is the bare electron mass, £l F c is the condensate volume, e F is the Fermi energy, and 

Q F is the volume of the Fermi sphere. One can imagine that the dispersionless plateau e(p) = fi 

given by Eq. (4) is slightly tilted counter-clockwise about fi and rounded off at the end points. 

If T <C Tf, it follows from Eqs. (1) and (5) that the effective mass M FC related to the FC is 

temperature dependent 

M*„ N(C\\ T, 

(7) 



M* FC N(0) T f 



M iVo(O) T ' 

where iVo(O) is the density of states of the noninteracting electron gas, and N(0) is the density 
of states at the Fermi level. We note that outside the FC region, the single-particle spectrum is 
not distinctly affected by temperature, being determined by the effective mass M£ given by Eq. 
(3), which is now evaluated at p < pi. Thus, we are led to the conclusion that systems with a 
FC must be characterized by two effective masses: M* FC related to the single particle spectrum 
of a low-energy scale and M£ related to the spectrum of a higher energy scale. The existence of 
these two effective masses can be observed as a break in the quasiparticle dispersion. This break 
is observed at temperatures T <^ Tf, and also when the superconducting state is superimposed 
on the FC state. In the former case, the occupation numbers over the area occupied by the 
fermion condensate are slightly disturbed by the pairing correlations such that the effective mass 
M FC becomes large but finite. We remark that at comparatively low temperatures, the FC and 
superconductivity go together because of the remarkable peculiarities of the FC phase transition. 
This transition is related to a spontaneous gauge symmetry breaking: the superconductivity order 



parameter n(p) = Jn(p)[l — n{p)\ has a nonzero value over the region occupied by the fermion 



condensate, while the gap A can vanish ||15 , |16 
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It is seen from Eq. (4) that at the FC phase transition point pf — > p i — > p Fj while the effective 
mass and the density of states tend to the infinity as follows from Eqs. (4) and (7). One can 
conclude that the beginning of the FC phase transition is related to the absolute growth of M FC . 
The onset of the charge density-wave instability in an electron system, which occurs as soon as the 
effective electron-electron interaction constant r s reaches its critical value r cdw , must be preceded 



by the unbounded growth of the effective mass [18]. For a simple electron liquid, the effective 



constant is proportional to the dimensionless average distance r s ~ ro/a_B between particles of 
the system in question, with r$ being the average distance and ag the Bohr radius. The physical 
reason for this growth is the contribution of the virtual charge density fluctuations to the effective 
mass. The excitation energy of these fluctuations becomes very small if r s ~ r cdw . Thus, a FC 
can occur when r s ~ r c d w . The standard Fermi liquid behavior can therefore be broken by strong 
charge fluctuations when the insulator regime is approached in a continuous fashion. We recall 
that the charge- density wave instability occurs in three-dimensional |19|] and two-dimensional (2D) 
electron liquids |2(J at a sufficiently high r s . As soon as r s reaches its critical value r FC < r cdw , 
the FC phase transition occurs. Thereafter, the condensate volume is proportional to r s — r F c and 



also Tf/sp ~ r s — r FC [jT^, 118] . In fact, the effective coupling constant r s increases with decreasing 
doping. It is assumed that both Tf and condensate volume VLpc build up with decreasing doping. 
The FC then serves as a stimulating source of new phase transitions lifting the degeneracy of the 
spectrum. The FC can produce, for instance, the spin density wave (SDW) phase transition or 
the antiferromagnetic one, thereby promoting a variety of the system properties. We note that the 
SDW phase transition, the antiferromagnetic transition, and the charge density one also depend on 
r s and occur at a sufficiently large value of r s even if the FC is absent. The superconducting phase 
transition is also aided by the FC. We analyze the situation where the superconductivity wins 
the competition with the other phase transitions up to a temperature T c . Above the temperature 
T* <C Tf the system under consideration is in its anomalous normal state, Eq. (7) is valid, and 
one can observe smooth non- dispersive segments of the spectra at the Fermi surface M . 



3 Superconductivity in the presence of FC 

We focus our attention on investigating the pseudogap that is formed above T c in underdoped 
(UD) high-temperature superconductors [4 — 8]. As we see in what follows, the existence of 
the pseudogap is closely allied with the presence of the FC characterized by a sufficiently high 
temperature Tf given by Eq. (6). Thus, the pseudogap is peculiar to UD samples, while optimally 
doped (OP) and overdoped (OD) samples may not exhibit this feature. We consider a 2D liquid 
on a simple square lattice that has a superconducting state with the d-wave symmetry of the order 
parameter k. We assume that the long-range component VJ r (q) of the particle-particle interaction 
Vp P (q) is repulsive and has the radius q\ r in the momentum space such that Pf/qw < 1- The 
short-range component V sr (q) is relatively large and attractive, with its radius Pf/Qst 1- In 
agreement with the d-symmetry requirements the low temperature gap A is then given by the 
expression 



21, 22, 23 



A(0) = 2k(4>)E((P) ~ Ax cos 20 = A^a; 2 - y 
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where E((f>) = ye 2 (0) + A 2 (0) and A x is the maximal gap. At finite temperatures, the equation 
for the gap can be written as 



p2n p 

A(p,0) = ~ J J ^ PP (p,0,Pi,0i)ft:(pi,0i)tanh 



2T Ait 2 



where p is the absolute value of the momentum and is the angle. It is also assumed that the 
FC arises near the van Hove singularities, leading to a large density of states at these points in 



accordance with Eq. (7). We note that the different FC areas overlap only slightly fL7f . A(</>) 
obeys the following equation that is determined by the chosen interaction V pp 

A(^ + 0) = -A(^ -</>). (9) 

It vanishes at tt/A and can therefore be expanded in the Taylor series around 7r/4, with p « pp\ 

A(p,0) = 9a -6 3 b + (10) 

where 9 = <f> — |. Hereafter, we consider solutions of Eq. (8) on the interval < 9 < n/A. 
We transform Eq. (8) by setting p « p F and separating the contribution l[ r coming from Vi r , 
with the contribution related to V sr denoted by I sr . At small angles, l[ T can be approximated 
in accordance with (10) by l\ T = 9 A + 9 3 B, with the parameters A and B independent of T if 
T <T* <^ Tf, because they are defined by the integral over the regions occupied by the FC. This 
theoretical observation is consistent with the experimental results showing that Ai is essentially 
T-independent at the temperatures T < T* |5]]. The coefficients of the expansion of I sr in powers 
of 9 depend on T. It is therefore more convenient to use the integral representation for I sr following 
from (8). We, thus, have 

A(0) = I sr + Iir = ~fjj Vsr(9, pi, 0i) tanh +9A + 9 3 B. (11) 

In Eq. (11), the variable p was omitted since p ~ pp. It is seen from this equation that the FC 
produces the free term 9 A + 9 3 B. In what follows, we show that at T > T no d e , the solution of Eq. 
(11) has the second node at 9 C (T) in the vicinity of the first node at tc/A. We also demonstrate 
that the temperature T noi } e has the meaning of the temperature T c at which the superconductivity 
vanishes. To show this, we simplify Eq. (11) to an algebraic equation. We have I sr ~ (Vq/T)9 
because tanh(£'/2T) ~ E/2T for E < T and T w T no ^ e , as is the case in the vicinity of the gap 
node at 9 = 0. The integration in Eq. (11) runs over a small area located at the gap node because 
of the small radius of V sr . Dividing both parts of Eq. (11) by k(9), we obtain 

E(9) = -(^-A l -9 2 B 1 )\9\, (12) 

where A\ and B\ are new constants and Vq ~ V, r (0) is a constant. Imposing the condition that 
Eq. (8) has the only solution A = when V sr = 0, we see that A\ is negative and B\ is positive. 
The factor in the brackets on the right-hand side of Eq. (12) changes its sign at some temperature 
T node w Vol A x] on the other hand, the excitation energy must be E{9) > 0. Therefore, we have 
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two possibilities p4| , p5| . The first follows from the assumption that A(9) = if 9 belongs to the 



interval VL n [0 < 6 < 9 C ]. In this case, for T > T node we must solve Eq. (8) with the condition 

A(0) = O, if : < 6? < 6» c ; T node <T. 

This resembles Eq. (4) with the parameter /z being equal to zero. The similarity is not coincidental, 
because we are searching for new solutions in both cases. Such solutions do exist because the points 
9 = and 9 = 9 C represent the branching points of the solutions. The second possibility can occur 
if the above solution does not lead to a minimum value of the free energy. Because the excitation 
energy must be positive for a stable state, the sign of A must be reversed at the point 9 = 9 C . 
Then the gap A(9) has the same sign within the interval Q n and changes its sign once more at the 
point 9 = 0, with A(9 C ) = A(0) = 0. Thus, we conclude that the gap A possesses new nodes at 
T > T no de HH) see Fig. 1. It can be seen from Eq. (12) that the angle 9 C is related to T > T no d e 



by 

TK (i™ (13) 

It follows from the above consideration and Eq. (12) that even below T no d e , the order parameter 
A cannot be approximated by a simple d-wave form; a more sophisticated expression must be 
used to fit the flattening of the gap A around the node. The following expression can be used for 
this purpose, 

A(0) = Ai[B cos(20) + (1 - B) cos(60)]. (14) 

Here < B < 1 in accordance with the experimental results and the term involving cos (60) 
is the next compatible with the d-symmetry of the gap. It also follows from Eq. (12) that the 
parameter B is a decreasing function of the temperature. At the temperatures T > T no ^ e , the 
value of (1 — B) is sufficiently large to produce new nodes of A given by Eq. (14). 

As an example of the solutions of Eqs. (8) and (11), we show, in Fig. 1 the gap A(0) 
calculated at three different temperatures 0.9T no d e , T no d e , and 1.2T noc fe. An important difference 
between curves (b) and (a) is the flattening of curve (b) at the nodes localized within the region 
Q n containing the interval —9 C < 9 < 9 C . As seen from Fig. 1, the flattening occurs as the result 
of the new nodes restricting the area fl n . It is also seen from Fig. 1 that the gap A is extremely 
small over the range Q n . 



It was recently shown in a number of papers (see, e.g., p(| £7|) that there exists an interplay 
between the magnetism and the superconductivity order parameters, leading to the damping of 
the magnetism order parameter below T c . Conversely, one can anticipate the damping of the 
superconductivity order parameter by magnetism. Thus, we conclude that the gap in the range 
fl n can be destroyed by strong antiferromagnetic correlations (or by spin density waves) existing 



in underdoped superconductors |^B, 2^]. It is believed that impurities can easily destroy A in the 
considered area. As a result, one is led to the conclusion that T c ~ T no d e , with the exact value of 
T c defined by the competition between the antiferromagnetic correlations (or spin density waves) 
and the superconducting correlations over the range fi n . 

We now consider the possibility for two quite different properties, the superconductivity and 
static spin density wave (SDW), to coexist. We start by briefly outlining the main features of the 
SDW p0| . A simple example is given by the linear SDW, with the net spin polarization P(r) 



P(r) = P ecos{Qx), (15) 
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radian 



Figure 1: The gap A as a function of <fi calculated at three different temperatures expressed in 
terms of T node ~ T c , while A is presented in terms of T*. Curve (a), solid line, shows the gap 
calculated at temperature 0.9T no d e - In curve (b), dashed line, the gap is given at T no d e - Note the 
important difference in curve (b) compared with curve (a) due to a flattening of the curve (b) 
over the region Q n . Calculated A(0) at 1.2T no d e is shown by curve (c), dotted line. The arrows 
indicate the two nodes restricting area fi n and emerged at T no d e - 



where Qx is the angle between the vectors Q and x. For convenience, the direction of the SDW 
is taken along the x axis, and e is the unit polarization vector, which in general can have any 
orientation with respect to Q. In contrast to the superconductivity, SDW can occupy only a part 
of the Fermi sphere with the volume 5S ~ pp5<f)5k, where 5<p is the Fermi surface angle and 5k is 
the "penetration depth" of the SDW into the Fermi sphere. At T = 0, the energy gain 5W due 
to the onset of SDW is given by 

SW ~ g 2 N(0)5(f>, (16) 



where g is the SDW gap determined by the formula f30| 
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PF<5k 

W) 



4 



exp( 



N(0) 7o 6<j> J 



where 70 is the coupling constant. As seen from Eq. (8), the variation of the gap within some area 
produces a variation of the gap over the entire occupied area with the same order of magnitude. 
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Therefore, elimination of A over a segment 5(f> requires the energy 8E\ ~ iV(0) A 2 (</>). We conclude 
that at T < T noi i e , the destruction of the gap on the interval 5(f> eliminates A over the entire region, 
because SEi is comparable with the gain 5E due to the superconducting state. A different situation 
occurs at the temperatures T > T no ^ when A is extremely small in Q n and the corresponding 
destruction energy satisfies inequality 5Ei -C 5E. Equations (16) and (17) are very similar to 
the corresponding BCS equations and this similarity also remains at finite temperatures |3(J 



Thus, the gain 8W and the gap g vary with the temperature similarly to the superconducting 
gain SE and the gap A. We also assume that the SDW transition temperature T n is sufficiently 
high, namely, T n > T c . We then come to the conclusion that 5Ei < 5W, and the region Q n is 
therefore occupied by the SDW at temperatures T > T node , resulting in the destruction of the 
superconductivity |24|, |25|| . We note that the Fermi surface angle 5(f> must be sufficiently large, 
because the gap g depends exponentially on 8(f) in accordance with Eq. (17). On the other hand, 
because we are dealing with SDW, we have 5(f>/ii ~ 1CU 2 |30fl . We thus conclude that a strong 
variation of the superconductivity characteristics may be observed in the vicinity of T node . 

It follows from the above considerations that A(#) can be destroyed only locally within the 
region Q n because of the different reasons. It also follows that T node is the temperature at which 
the superconductivity vanishes, that is, T c ~ T no d e - As to the gap at T > T c , or more precisely, 
the pseudogap, it persists outside the Q n region. In accordance with [|, fjj, we see that the 
superconducting gap A(#) smoothly transforms into the pseudogap at T > T c . We can therefore 
expect a dramatic reduction in the difference between the free energy of the normal and the 
superconducting state at T = T c (the so-called condensation energy, which we consider in some 
detail in the next section). It can then be concluded that the temperature T* has the physical 
meaning of the BCS transition temperature between the state with the order parameter k / 



and the normal state. Because T c ~ Vq/A\, we find from Eq. (13) that 6 C ~ (T — T c )/T c . This 
result is in harmony with our calculations of the function 6 C ([T — T c ]/T c ) plotted in Fig. 2. Thus, 
we conclude that the pseudogap "dies out" in UD samples as the temperature T* is approached. 
Quite naturally, one has to recognize that Ai must scale with T* . 

A few remarks are in order at this point. On the basis of the previous consideration, we conclude 
that the BCS approach is fruitful in considering OD, OP, and UD samples in the weak coupling 
regime. With more underdoping, the antiferromagnetic correlations become stronger, breaking 
down the gap over the range fl n at lower temperatures. Thus, one observes the decrease of T c 
with the decrease of doping. On the other hand, the condensate volume Qfc becomes larger with 
the decrease of doping, leading to increase of the gap Ai which is proportional to the volume and 
interaction V pp [0. Consequently, the temperature T* becomes higher with decreasing doping. 
All these results are in agreement with the experimental findings @, Q. A peak was observed at 
41 meV~ 2Ai in inelastic neutron scattering from single crystals of the OD, OP, and UD samples 
YBa 2 Cu 3 6+:1 . and Bi 2 Sr2CaCu208+,5 at temperatures below T c , while a broad maximum above T c 
exists in underdoped samples only |5T, [!2|. The explanation of this peak given in |33[ was based 



on the ideas of the BCS theory. From the above discussion, it appears that the same explanation 
holds for the broad maximum in underdoped samples above T c because the physics of the process 
is essentially the same. 
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Figure 2: Calculated angle 9 C , pulling apart the two nodes, as a function of (T — T c )/T c . 

4 Condensation energy 

We now consider the energy gain or condensation energy E con d liberated when the system in 
question undergoes the superconducting phase transition involved in the FC phase transition. We 
set T = for simplicity. The energy E con d can be schematically broken into two parts related 



to the kinetic and the potential energy. The condensation energy was considered in [p4 |, where 
it was argued that the main contribution to the condensation energy comes from the kinetic 
energy, i.e., the superconducting phase transition of high-temperature superconductors is kinetic 
energy driven. Here, we give a possible interpretation of the situation. It is known |35[] that in 
the superconducting phase transition, the positive contribution comes from the potential energy, 
while the gain in the kinetic energy is negative. In the other words, the superconducting phase 
transition is driven by the gain in the potential energy. This result is rather obvious because the 
ground state energy E gs is given by 

E gs [K{p)) = E[n{p)) + E sc [k{p)), (18) 



with the occupation numbers n(p) determined by /t(p) = yn(p)[l — n(p)]. The second term 
E sc [k(p)] on the right-hand side of Eq. (18) is defined by the superconducting contribution, which 
in the simplest case is of the form 



E sc [k(p)] = 92 J ^pp(pi,p 2 )/«(pi)/«(p2) 



dp 1 dp 2 



The first term i?[ra(p)] can be taken as 

2 dp , g x 



E[n{p)] = I 2M n{p) 4^ + 9 iJ ^ fa' P*)»(PiMPa 



dpi<ip2 



(19) 



(20) 
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with the second integral playing the role of the exchange-correlation contribution to the ground 
state energy. If the effective mass M* L given by Eq. (3) is positive and finite, i£[n(p)] reaches its 
minimum at n(p) = n F (p) and increases with the deviation of n(p) from the Fermi distribution, 
as it occurs in the presence of superconducting correlations. Thus, the standard situation is 



that the superconducting phase transition is driven by a decrease of the potential energy [35 
The situation can be different if the system undergoes the FC phase transition. To see this we 
temporarily assume that cfe — ► and rewrite Eq. (20) as 

£7[n(p)] = / e(p)n(p)0 - f / V( Pl , p 2 )n( Pl )n(p 2 )^^, (21) 

with the single particle energy 

6E[n(p)] 

£(p) = -sMpf (22) 

The energy £7[n(p)] can be lowered by alteration of n(p) if Eq. (4) has solutions. As the result, 
we can write the inequality JIT 



E cond = E N -E FC > J (e(p) - fi) Sn(p)^ > 0, (23) 

with En being the energy of system in its normal state, Epc the energy with FC, and the integral 
taken over the region occupied by FC. The chemical potential \x preserves the conservation of 
the particle number under the variation 5n(p). We assume that the kinetic energy is given by 
the first term on the right-hand side of Eq. (21). It then follows from Eq. (23) that the kinetic 
energy can be lowered, and this lowering is driven by the FC phase transition. It is instructive to 
illustrate this by a simple example. We take V(pi,p2) = gi<5(pi — P2), then E con d given by Eq. 
(23) becomes 

E cond = J MPK(P) - e(p)n(p)) ^ + f / (n 2 (p) - n F (p)) ^ (24) 

with £o(p) being the single particle energy of the normal ground state. It is easily verified that 
the second term on the right-hand side of Eq. (24), which is related to the potential energy gain, 
is negative. This term can be written as 

y / ( n2 (p) ~ n Up)) = % J ~ n F(p)) {n{p) + n F (p)) 
Observing that 

(n(p) - n F {p)) ^ = 
because of the particle number conservation and taking into account that 

{n{p) +n F (p))\ p < PF > (n(p) + n F (p))\ PF < pi 

we arrive at the conclusion. The first term is positive because of inequality (23). Thus, we are led 
to the conclusion that the FC phase transition can be considered as driven by the kinetic energy. 
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We now let the coupling constant g 2 be small, then the gap A is proportional to g 2 pT| . The 
optimum values of the occupation numbers given by Eq. (4) are disturbed, leading to an increase 
of the energy i?[n(p)]. The positive gain in the potential energy given by Eq. (19) is driving the 
formation of the superconducting ground state. Because the coupling constant g 2 is sufficiently 
small, the structure of the system ground state is defined by the FC, and the superconducting 
state is a "shadow" of the FC under these conditions ||15|| . Then, the main contribution to E con d 



comes from the FC phase transition, and the complex transition (FC plus superconductivity) is 
kinetic energy driven [BBJ. On the other hand, in the case where FC is weak compared to the 
superconductivity (or is absent), we are dealing with a pure superconducting phase transition, 
which is obviously potential energy driven. 



5 Quasiparticle dispersion and lineshape 

We now discuss the origin of two effective masses M£ and Mp C occurring in the superconducting 
state and leading to a nontrivial quasiparticle dispersion and a change of the quasiparticle velocity. 
As we see in what follows, our results are in a reasonably good agreement with the experimentally 



deduced data flTO, |L El. For simplicity we set T = 0. Varying E gs given by Eq. (18) with respect 



to a p , we find 

E gA a P } 
5a p 



(e(p) - /i) tanh(2a p ) + A(p) = 0, (25) 



with n(p) = cos (« p ), n(p) = sin(a p ) cos(a p ), and e(p) defined by Eq. (22). As g 2 — > 0, we have 
that A(p) — > 0, and Eq. (25) becomes 

(e(p) - /i) tanh(2a p ) = 0. (26) 

Equation (26) requires that 

e (p) -fi=0, if tanh(2a p ) ^ (0 < n(p) < 1), (27) 



which leads to the FC solutions defined by Eq. (4) |16], Pq| . As soon as the coupling constant 
g 2 becomes finite but small, such that g 2 /g± <C 1, the plateau e(p) — /x = is slightly tilted and 
rounded off at the end points. This implies that 

£ (p)-/i~A 1 , (28) 

which allows us to estimate the effective mass as 



M* FC T f 



M Ai 



(29) 



Outside the condensate area, the quasiparticle dispersion is determined by the effective mass M£ 
given by Eq. (3). We note that calculations in the context of a simple model support the above 
consideration |15|| . In that case, putting V(pi,p2) = 5(pi,p2) and V^, p (pi,p2) = 5(pi,p2) in Eqs. 
(19) and (20) and carrying out direct calculations, we obtain at T = 

E = e( Pf ) - e( Pi ) ~ {Pf ~/: )PF * 2A X . (30) 

IV1 FC 
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On the other hand, at T > T c , taking into account that n(pi) ~ 1 and n(pf) ~ 0, we obtain from 
Eq. (5) with the same accuracy, 

E * (P/-^ * 2T. (31) 

Equations (30) and (31) allow us to estimate the effective mass M FC related to the region occupied 
by the FC at temperatures T ^ Tf. Outside the region, the effective mass is M£. When Eqs. 
(28) and (29) are compared with Eqs. (5) and (7), it is apparent that the gap Ai plays the role 
of the effective temperature that defines the slope of the plateau. On the other hand, at T = T c 
in OD or OP samples, the gap vanishes and Eqs. (5) and (31) define the quasiparticle dispersion 
and the effective mass. Taking into account that A x ~ T c , we are led to the conclusion that Eqs. 
(28) and (29) derived at T = match Eqs. (5) and (7) at T c . Thus, Eqs. (28) and (29) are 
approximately valid over the range < T < T c . It follows from Eq. (30) that at T < T c , the 
quasiparticle dispersion can be presented with two straight lines characterized by the respective 
effective masses M FC and M£ and intersecting near the binding energy Eq ~ 2Ai. Equation (31) 
implies above T c , the lines intersect near the binding energy ~ 2T. The break separating the faster 
dispersing high-energy part related to M£ from the slower dispersing low-energy part defined by 
M FC is likely to be enhanced in UD samples at least because of the rise of the temperature Tf, 
which grows with the decrease of doping. We recall that in accordance with our assumption, 
the condensate volume Qfc and Tf are growing with underdoping, see Eq. (6) and Sec. III. 
It was also suggested that the FC arises near the van Hove singularities, while the FC different 
areas overlap only slightly. Therefore, as one moves from (0,0) towards (vr,0) the ratio M FC /M^ 
grows in magnitude, developing into the distinct break. In fact, assuming that the temperature 
Tf depends on the angle along the Fermi surface and taking Eq. (29) into account, one can 
arrive at the same conclusion. The dispersions above T c exhibit the same structure except that 
the effective mass M FC is governed by Eq. (31) rather than (30) and both the dispersion and the 
break are partly "covered" by the quasiparticle width. Thus, one concludes that there also exists 
a new energy scale at T <C Tf defined by E and intimately related to Tf [[36| . 



We turn to the quasiparticle excitations with the energy E{4>) = ye 2 (0) + A 2 (0). At temper- 
atures T < T c , they are typical excitations of the superconducting state. We now qualitatively 
analyze the processes contributing to the width 7. Within the limits of the analysis, we can take 
A ~ 0, which corresponds to considering excitations at the node. Our treatment is then valid for 
both T <T C and T c < T. For definiteness, we consider the decay of a particle with the momentum 
p > pf- Then j(p,u) is given by the imaginary part of the diagram shown in Fig. 3a, where 
the wiggly lines stand for the effective interaction. Because of the unitarity, diagram 3b which 



represents the real events can be used to calculate the width [37 



7(p,w) = 2tt 



V(q) 



i(q, -uJ pq ) 



n(k)(l-n(k + q))i(w M + w fc! )— , (32) 



with e(q, —uj pq ) being the complex dielectric constant and V(q)/e the effective interaction. Here, 
q and u>k q = e(k + q) — e(k) are the transferred momentum and energy, respectively, and u pq = 
to — e(p — q) is the decrease in the quasiparticle energy as the result of the rescattering processes: 
the quasiparticle with the energy ui decays into a quasihole e(k) and two quasiparticles e(p — q) 
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p-q 



P 



P 



p-q 





k 



k 



a 



Figure 3: Diagram 3a depicts a process contributing to the imaginary part. Diagram 3b shows a 
real process contributing to the imaginary part, observe that quasiparticles p — q, k + q and k 
are on the mass shell. 



and e(k + q). The transferred momentum q must satisfy the condition 



V > P ~~ Q > Pf- 



(33) 



Equation (32) gives the width as a function of p and uj; the width of a quasiparticle with the 
energy e(p) is given by j(p, uj — e(p)). 



Estimating the width in Eq. (32) with the constraint (33) and uj pq ~ T, we find that 



7 (p,u; = e{p)) ~ {MlfT 2 



(34) 



for normal Fermi liquids. In the case of the FC one could estimate 7 ~ 1/T upon using Eqs. (9) 
and (34). This estimate were correct if the dielectric constant is small, but e 
result, for the FC we have 



M* FC . As the 



7(p,w = e{p)) 



(M 



FC j 



\3rp2 



FC) 



e F 



(35) 



where Ef is the Fermi energy | 38fl . Calculating 7Q9, uj) as a function of p at constant uj, we obtain 
the same result for the width given by Eq. (35) when uj = e(p). The calculated function can 
be fitted with a simple Lorentzian form, because quasiparticles and quasiholes involved in the 
process are also located in the vicinity of the Fermi level provided uj — Ef ~ T. It then follows 
from Eq. (35) that the well-defined excitations exist at the Fermi surface even in the normal 
state |3"5fl. This result is in line with the experimental findings determined from the scans at a 



constant binding energy (momentum distribution curves or MDCs) [|10], |39|. On the other hand, 
considering j(p,u) as a function of uj at constant p, we can check that the quasiparticles and 
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quasiholes contributing to the function can have the energy of the same order of the magnitude. 
For uj — Ep ~ T, the function is of the same Lorentzian form, otherwise the shape of the function is 
disturbed at high uj by high-energy excitations. In that case the special form of the quasiparticle 
dispersion characterized by the two effective masses must be taken into account. As the result, 
the lineshape of the quasiparticle peak as a function of the binding energy possesses a complex 
peak-dip-hump structure || |9|, [5J directly defined by the existence of the effective masses M FC 
and M£. Our consideration shows that it is the spectral peak obtained from MDCs that provides 
important information on the existence of well-defined excitations at the Fermi level and their 



width p6[ . The detailed numerical results will be presented elsewhere. 

At T > T c , the gap is absent in OD or OP samples, and the width 7 of excitations close 
to the Fermi surface is given by Eq. (35). For UD samples, A(#) =0 in the range Q n and we 
have normal quasiparticle excitations with the width 7. Outside the range Q n , the Fermi level is 
occupied by the BCS-type excitations with finite excitation energy given by the gap A(#). Both 
types of excitations have widths of the same order of magnitude. We now estimate 7. For the 
entire Fermi level occupied by the normal state, the width is equal to 7 ~ N(0) 3 T 2 /(3 2 , with 
the density of states N(0) ~ 1/T and the dielectric constant (3 ~ N(0). Thus, 7 ~ T |15|]. In 
our case, however only a part of the Fermi level within Q n belongs to the normal excitations. 
Therefore, the number of states allowed for quasiparticles and for quasiholes is proportional to 
9 C , the factor T 2 is therefore replaced by T 2 9 2 . Taking these factors into account, we obtain 
7 ~ 9 2 T ~ T(T — T c )/T c ~ (T — T c ), because only small angles are considered. Here, we have 
omitted the small contribution coming from the BCS-type excitations. That is why the width 
7 vanishes at T = T c . Thus, the foregoing analysis shows that in UD samples at T > T c , the 
superconducting gap smoothly transforms into the pseudogap. The excitations of the gapped area 
of the Fermi surface have the same width 7 ~ (T — T c ) and the region occupied by the pseudogap 
is shrinking with increasing temperature. These results are in good qualitative agreement with 
the experimental facts [4 — 7] . 



6 Concluding remarks 

We have discussed the model of a strongly correlated electron liquid based on the FC phase 
transition and extended it to high-temperature superconductors. The FC transition plays the role 
of a boundary separating the region of a strongly interacting electron liquid from the region of a 
strongly correlated electron liquid. On the basis of the BCS theory ideas we have also considered 
the superconductivity with the d-wave symmetry of the order parameter in the presence of the 
FC. We can conclude that the BCS-type approach is fruitful for OD, OP, and UD samples. We 
have shown that in UD samples, the gap becomes flatter near the nodes at temperatures T < T c , 
and the superconducting gap smoothly transforms into a pseudogap above T c . The pseudogap 
occupies only a part of the Fermi surface, which eventually shrinks with increasing temperature, 
vanishing at T = T*, and the maximum gap A x scales with the temperature T*. We have also 
shown that the general dependence of T c , T*, and Ai on the underdoping level fits naturally into 
the considered model. At temperatures T* > T > T c , the single-particle excitations of the gapped 
area of the Fermi surface have the width 7 ~ (T — T c ). The quasiparticle dispersion in systems 
with FC can be represented by two straight lines characterized by the respective effective masses 
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Mp C and M£. At T < T c , these lines intersect near the point Eq ~ 2Ai, while above T c , we 
have E'o ~ IT . It is argued that this strong change of the quasiparticle dispersion at Eq can be 
enhanced in UD samples because of strengthening the FC influence. The single-particle excitations 
and their width 7 are also studied. We have shown that well-defined excitations with 7 ~ T exist 
at the Fermi level even in the normal state. This result is in line with the experimental findings 
determined from the scans at a constant binding energy, or MDCs. We have also treated the FC 
phase transition in the presence of the superconductivity and shown that this phase transition can 
be considered as kinetic energy driven. Thus, without any adjustable parameters, a number of the 
fundamental problems of strongly correlated systems are naturally explained within the proposed 
model. 

This research was supported in part by the Russian Foundation for Basic Research under Grant 
No. 98-02-16170. 
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